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We present a comprehensive examination of the X-ray variability of the narrow line Seyfert 1 (NLSl) galaxy 
NGC 4051. We combine over 6.5 years of frequent monitoring observations by RXTE with a > 100 ks continuous 
observation by XMM-Newton and so present a powerspectral density (PSD) covering an unprecedent frequency 
range of over 6.5 decades from < 10"* to > 10"^ Hz. The combined RXTE and XMM-Newton PSD is a very 
good match to the PSD of the galactic black hole binary system (GBH) Cyg X-1 when in a 'high', rather than 
'low', state providing the first definite confirmation of an AGN in a 'high' state. 

We find a break in the PSD at a frequency ub = 8I3 x 10~'*Hz. If 1/3 scales linearly with mass then, assuming a 
black hole (BH) mass of WMq for Cyg X-1, we imply a BH mass of Stl x 10^ Mq in NGC 4051, which is consistent 
with the recently reported reverberation value of 5^3 x W^Mq. Hence NGC 4051 is emitting at ~ 30% LEdd. 

We note that the higher energy photons lag the lower energy ones and that the lag is greater for variations of 
longer Fourier period and increases with the energy separation of the bands. Variations in difi'erent wavebands 
are very coherent at long Fourier periods but the coherence decreases at shorter periods and as the energy 
separation between bands increases. This behaviour is similar to that of GBHs and suggests a radial distribution 
of frequencies and photon energies with higher energies and higher frequencies being associated with smaller radii. 

It is not possible to fit all AGN to the same linear scaling of break timescale with black hole mass. Broad line 
AGN are consistent with a linear scaling of break timescale with mass from Cyg X-1 in its low state but NLSl 
galaxies scale better with Cyg X-1 in its high state. We suggest that the relationship between black hole mass 
and break timescale is a function of another underlying parameter which may be accretion rate. 



1. INTRODUCTION 

Early observations with EXOSAT, coupled 
with some archival data, indicated a similarity 
between the powerspectral densities (PSDs) of 
AGN and galactic black hole X-ray binary sys- 
tems (GBHs). Both types of PSD were steep 
at high frequencies and flattened, below a break 
or 'knee' frequency |1(J|. Early indications were 
that the break timescale scaled with black hole 
mass but the long timescale (weeks - years) X-ray 
lightcurves of AGN were not of adequate quality 
to properly determine the break timescales. How- 
ever since the launch of the Rossi X-ray Timing 
Explorer {RXTE; JW) in December 1995 we have 
been able to obtain excellent long timescale AGN 



lightcurves, and PSDs, and evidence is building 
for a correlation between black hole mass and 
break timescale (eg |11I3I8| ). 

GBHs, however, are found in at least two 
states: the low flux, hard spectrum ('low') state 
and the high flux, soft spectrum ('high') state (eg 
I9II3I2) ). The PSDs of these two states are differ- 
ent. At high frequencies both PSDs have slopes of 
^ —2.5. Below a break frequency (vb 3Hz for 
the low state and ~ 15Hz for the high state) the 
PSDs both flatten to a slope of -1. In the high 
state there is no other observable slope change 
but in the low state the PSD flattens further 
(below ly ^ 0.3Hz) to a slope of zero. It had 
been thought, based on the similarity of their en- 
ergy spectra, that AGN were the analogues of 
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Figure 1. RXTE Long Term 2-10 keV lightcurve of 
NGC 4051. Each data point represents an observation 
of ~ 1 ks. 



GBHs in their low state. However, by combining 
RXTE and XMM- Newton observations, we show 
that NGC4051 is the analogue of a high state 
GBH. The PSD we present of NGC4051 is the 
best PSD yet produced for an AGN. In this paper 
we consider the implications of our result for the 
understanding of the similarity between AGN and 
GBHs. We also examine the variation, with en- 
ergy, of the high frequency PSD slope and the 
coherence of the X-ray emission, and draw con- 
clusions regarding the geometry of the emitting 
region and the origin of theblack holemass/break 
timescale relationship. 

2. OBSERVATIONS 

In Fig. n we present all of the short (^1 ks 
duration) monitoring observations of NGC 4051 
which have been carried out with the proportional 
counter array (PGA, on RXTE from 1996 

until 2002 inclusive. In Fig.|21we present the full- 
band XMM-Newton lightcurve obtained from a 
continuous observation over the 16th and 17th 
May 2001. The lightcurve is the sum of data from 
the PN and MOS detectors. We refer readers to 
|12| for details of the analysis of the RXTE and 
XMM-Newton data. 
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Figure 2. XMM-Newton background-subtracted 

lightcurve of NGC 4051 in the 0.1-10 keV energy band, 
with 5s time bins. 



3. POWERSPECTRAL ANALYSIS 

3.1. Combined RXTE and XMM- 
Newton PSD 

Powerspectral analysis of the continuous 
XMM-Newton lightcurves is straighforward but 
analysis of the RXTE lightcurves is more diffi- 
cult as they consist of many short observations 
with a variety of sampling patterns in order to 
cover different Fourier periods. Since early 2000 
we have had observations approximately once ev- 
ery 2 days, and we have had a period (mid 2000) 
of 2 months where we observed every 6 hours. 
However prior to 2000 the observations were less 
frequent. To cope with the distortions which the 
sampling pattern introduces into the PSD, we em- 
ploy a simulation-based modelling technique |21j . 
We are able to incorporate the XMM-Newton ob- 
servations into our modelling. We again refer to 
[E] for details. Below (Sect ion ESI 

we show that 

vb is independent of energy and so here we use 
the 0.1-2 keV XMM-Newton data which shows 
the break most clearly. 

A simple powerlaw, although a reasonable fit 
to the low frequency RXTE PSD on its own, is a 
very poor fit to the combined RXTE and XMM- 
Newton PSD. The best fit (FigEJ is to a bending 
powerlaw ll2j where the PSD changes smoothly 
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Figure 3. Combined RXTE and XMM-Newton PSD. 
The underlying, undistorted, model is given by the smooth 
dashed line. The observed data is given by the continu- 
ous jerky line and the model, distorted by the sampling 
pattern, is given by the points with errorbars. 



from a slope of ~ —1.1 below, to a slope of 
aH = -2.8 above vb = 8^3 x lO^'^Hz. Note 
particularly the lack of a second break, about a 
decade below i^b , which is characteristic of a low 
state GBH. 

In order to provide a good template for compar- 
ison we present, in FigQ] the PSD of the GBH 
Cyg X-1 in its high state. It is fitted very well 
by the same bending powcrlaw model which fits 
Fig 13 We therefore conclude that NGC 4051 is 
the first AGN confirmed to be in a 'high' state. 

In Cyg X-1 (FigEJ vb = 22± IHz. If vb scales 
linearly with mass then, assuming a black hole 
mass of IOMq for Cyg X-1 0, we imply a black 
hole mass of 3^1 x IO^Mq in NGC 4051, which is 
consistent with the recently reported reverbera- 
tion value of 5^^ x IO^Mq [Hj. Hence NGC 4051 
is emitting at ~ 30% LEdd- 

3.2. Variation of high frequency XMM- 
Newton PSD with energy 

We made PSDs of NGC 4051 in a number 
of different XMM-Newton wavebands. All PSDs 
showed a break at i^b ^ 10~^ Hz and had similar 
slopes below i/b- As the PSD slope below i^b in 
Cyg X-1 in the high state is also independent of 
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Figure 4. PSD of Cyg X-1 in the high state in the 2- 
13 keV range. The dashed line is a powerlaw with an 
exponential cut-off. The solid line is a bending powerlaw, 
which is a marginally better fit. Note the very strong 
similarity to the PSD of NGC 451 (Fig. 



energy T we fixed at -1.1, the best fit value 
from Fig 13 Further fitting confirmed the invari- 
ance of vb with energy but showed a systematic 
flattening of an with increasing energy (Figl^l). 
Similar behaviour is seen in MCG-6-30-15 |23) . 

Flattening of the PSD with increasing energy 
at high frequency is seen in Cyg X-1 in the low 
state 13^. In the generally accepted Compton 
scattering scenario for the production of X-rays 
the flattening is generally interpreted as arising 
from variations in the scattering corona as the 
increased scattering necessary to raise photons to 
higher energies would wash out variations from 
the seed photons. 

4. THE RMS-FLUX RELATIONSHIP 

Uttley and M^Hardy 20 showed that the am- 
plitude of absolute rms variability in small seg- 
ments of the X-ray lightcurves of GBHs scaled 
remarkably linearly with the mean flux of the seg- 
ments. They also showed that the variability am- 
plitude in RXTE lightcurves of Seyfert galaxies 
increased with X-ray flux, but could not tell if the 
relationship was really linear. Subsequently, bet- 
ter quality lightcurves have provided stronger ev- 
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Figure 5. PSD of NGC 4051 as measured by XMM- 
Newton in the 0.1-2 keV energy band. The soUd Hne is 
the best fit with fixed at -1.1. The Poisson noise level 
has been subtracted from the PSD. The dashed line is the 
best 2-10 keV fit with normalisation adjusted so that the 
PSDs agree at low frequencies. 



idence of a linear rms-flux relation in AGN |4l2;^j . 

We measured the mean flux and the noise- 
subtracted variance of 39 separate segments, of 
duration 2560-s, of the 5-s binned 0.1-10 keV 
XMM-Newton NGC 4051 lightcurve. We then 
averaged the variance in 4 flux bins and took the 
square root to determine the rms (Fig. EJ. We 
note a strong linear rms-flux relationship, similar 
to that seen in GBHs, indicating a strong simi- 
larity in the process driving the variability. 

5. TIME LAGS AND COHERENCE 

Measurements of 'time lags' and 'coherence' be- 
tween different wavebands provide a strong di- 
agnostic of the geometry of the emission region 
(eg |22l(jj and Section l?rT|l . Here, for each wave- 
band, we split the lightcurve into its Fourier com- 
ponents and measure the lag between variations 
in different bands, as a function of Fourier period. 
The results are shown in Fig[7| The hard band 
always lags the soft band and the lag increases 
with increasing Fourier period and with increas- 
ing energy separation of the bands (cf IT5). 

The coherence between two lightcurves, at a 
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Figure 6. rms-flux relationship for NGC 4051 derived 
from the XMM-Newton observations in the 0.1-10 keV 
band. 



certain frequency, may be interpreted as the cor- 
relation coefficient between the Fourier compo- 
nents of the two lightcurves at that frequency 
Pl^ . For NGC4051 we note (FigEJ that, even 
taking account of the artificial decrease in co- 
herence introduced by numerical approximations, 
the coherence decreases with increasing band sep- 
aration and with increasing Fourier period. 

6. DISCUSSION 

6.1. Geometry of the Emission Region 

The linear rms-flux relationship tells us that 
the amplitude of short timescale variations re- 
sponds to the amplitude of longer timescale vari- 
ations. This response is explained naturally 
in the context of the model of Lyubarskii [7] 
where a characteristic timescale is associated with 
each radius of the accretion disc, with shorter 
timescales arising closer to the black hole. Vari- 
ations propogate inwards and modify the ampli- 
tude of shorter timescale variations. If we also 
associate different photon energies with different 
locations, with the bulk of the higher energies be- 
ing produced closer in, we can explain the obser- 
vations of lags and coherence mentioned earlier. 

Churazov et al. PP, with enhancements by 
Kotov et al. have built a qualitative model 




Fi gur6 7. Time lag vs. Fourier period for the cross 
spectrum of the soft energy band (0.1 — 0.5 keV) vs. the 
0.5 — 2, 2 — 5, and 5 — 10 keV bands (filled, open circles 
and filled squares, respectively). The hard bands are all 
delayed with respect to the soft (0.1-0.5 keV) band. The 
solid, dashed, and dot-dashed lines show the best-fitting 
power law models, assuming a slope of 1. Note that the 
open and filled circles have been displaced slightly in pe- 
riod so that errorbars do not overlap. 



Figure 8. Coherence function vs. Fourier period for the 
soft energy band (0.1 - 0.5 keV) vs. the 0.5 - 2 and 2-10 
keV bands (filled circles and open squares, respectively). 
The two open squares at the longest Fourier period have 
been displaced in period to avoid overlapping of errorbars. 
The dashed and dot-dashed lines show, for the 0.1 — 0.5 
keV vs. 0.5 — 2 and 2 — 10 keV bands respectively, the 95% 
confidence limit for spurious lack of coherence introduced 
by the approximation in the numerical calculation. 



along these lines. In the Kotov model variations 
propogate inwards until they hit the X-ray emit- 
ting region, whose emission they modulate. As a 
larger fraction of the lower energy emission is pro- 
duced further out in the disc, where characteristic 
timescales are longer, we expect larger lags with 
increasing energy separation of the bands and at 
longer Fourier periods. Regarding the coherence 
then, if the highest energies are preferentially pro- 
duced at the innermost radii, those energies see 
the largest spectrum of variations, including those 
at the shortest periods. However if the perturba- 
tions only propagate inwards then the lower ener- 
gies do not see the shortest period variations and 
the coherence then decreases as we go to shorter 
periods and greater separation of the bands. 



6.2. AGN black hole masses and break 
timescales 

We have gathered data from the literature (eg 
|8l21ll5l23l3^ ) and have plotted the timescale at 
which the PSD breaks from a slope of ~ —1 to a 
slope of ^ —2.5, against the best available value 
of the black hole mass. This break might be the 
only break in a high state PSD or it might be 
the higher frequency break in a low state PSD. 
For no AGN but NGC 4051 can we be sure which 
break it is. We also plot the break for Cyg X-1 in 
its high and low states. We plot broad line AGN 
with filled circles and NLSls with open circles. 
The best fit line to all of the AGN together (not 
shown) does not go through either state of Cyg 
X-1. 

The data are not sufficient to perform mean- 
ingful fits to subsets of the data but we can plot 
lines of linear scaling through both the high and 
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Figure 9. PSD break timescales vs black hole mass. 
NLSls are shown as open circles and broad line Seyferts 
are shown as filled circles. The high (H) and low (L) states 
of Cyg X-1 are also plotted and lines (dot-dash and long- 
dash respectively) of slope 1.0 are drawn through those 
points. The solid arrow labelled with fn indicates the way 
that the break timescale / mass line may move with increas- 
ing accretion rate. 



low states of Cyg X-1. We note that the broad 
Une AGN are consistent with hnear scaHng from 
Cyg X-1 in its low state but the NLSls scale bet- 
ter to the high state. We therefore suggest that 
some parameter varies between the high and low 
states in such a way as to move the line of the 
scaling relationship to shorter timescales for the 
NLSls/high state systems. The break timescale 
in NGC4051 (1250s) corresponds approximately 
to the viscous or thermal timescale at ~few Rg 
(cf mi). In GBHs the inner edge of the accretion 
disc moves inwards as the accretion rate increases 
(eg If the break timescale is therefore asso- 
ciated with the inner edge of the accretion disc, 
the parameter which moves the position of the 
mass-timescale scaling relationship might well be 
the accretion rate, m. 
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